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Abstract 

We extend Hardy's uncertainty principle for a square integrable 
function ip and its Fourier transform to the n-dimensional case using a 
symplectic diagonalization. We use this extension to show that Hardy's 
uncertainty principle is equivalent to a statement on the Wigner dis- 
tribution Wip of tp. We give a geometric interpretation of our results 
in terms of the notion of symplectic capacity of an ellipsoid. Further- 
more, we show that Hardy's uncertainty principle is valid for a general 
Lagrangian frame of the phase space. Finally, we discuss an extension 
of Hardy's theorem for the Wigner distribution for exponentials with 
convex exponents. 



1 Introduction 

A folk metatheorem is that a function tp and its Fourier transform Ftp 
cannot be simultaneously sharply localized. An obvious manifestation of 

*This author has been supported by the European Union EUCETIFA grant MEXT- 
CT-2004-517154. 
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this "principle" is when ip is of compact support: in this case the Fourier 
transform Ftp can be extended into an entire function, and is hence never of 
compact support. A less trivial way to express this kind of trade-off between 
ip and Ftp was discovered in 1933 by G.H. Hardy [15]. Hardy showed, using 
methods from complex analysis (the Phragmen-Lindelof principle), that if 
ip £ L 2 (M) and its Fourier transform 



oo 



Ftp(p) = 1 / e n px ^( x )dx 
satisfy, for \x\ + \p\ — > oo, estimates 

ij){x) = 0{e~m x2 ) , F%l)(p) = 0(e'^h p2 ) (1) 
with a, b > 0, then the following holds true: 

(1) If ab>l then tp = 0; 

(2) If ab = 1 we have ip{x) = Ce~2K x2 for some complex constant C ; 

(3) If ab < 1 we have ip(x) = Q(x)e~2h x where Q is a polynomial func- 
tion. 

Recently, researchers in harmonic analysis and time-frequency analysis 
have formulated variants of Hardy's theorem for phase space representations 
(time-frequency representations) such as the Wigner distribution, see 
[TOl Illj . The results in \10\ [TT] are deduced from Hardy's theorem for a 
carefully chosen square-integrable function and its Fourier transform. In 
[2] a multidimensional extension of Hardy's theorem is presented, which 
are based on an extension of the Phragmen-Lindelof principle to several 
complex variables. The results of [2] have in a sense the same flavor as 
our statements, but they are of a completely different nature. Actually, we 
only invoke real variable methods in our proof of the n-dimensional Hardy 
theorem. 

The principal aim of this paper is to reformulate Hardy's theorem in 
terms of phase-space objects. We will actually give a non-trivial restatement 
of Hardy's theorem for functions ip G L 2 (M. n ) satisfying estimates 

ip(x)=0(e-m xTAx ) , Fi,{p) = 0{e-^h xTBx ) (2) 

where A, B are positive-definite symmetric matrices, and show that the es- 
timates 

iP(x) =0{e-^h Ax2 ) , Fi>(p) = 0{e~^h Bx2 ) (3) 
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are equivalent to a single estimate 

Wip(x,p) = 0{e~Ti {xTAx+pTBp) ) for \x\ + \p\ -> oo 

for the Wigner transform of ip. This theorem provides a positive answer 
to a question raised by Grochenig in [11] on the equivalence of uncertainty 
principles for a function and its Fourier transform and uncertainty prin- 
ciples for the Wigner distribution (or more generally, for any phase space 
representation) . 

We will see that the geometric interpretation of the conditions on the 
matrices A, B is that the symplectic capacity of the "Wigner ellipsoid" 

W : x T Ax + p T Bp < h 

is at least ^h, the half of the quantum of action. This property is related to 
the fact that the notion of symplectic capacity is a natural tool for expressing 
the uncertainty principle of quantum mechanics in a symplectically covariant 
and intrinsic form as discussed in de Gosson [5j EJ [7] ; also see de Gosson 
and Luef [3 [9]); it turns out that, more generally, if a function ip £ L 2 (M™) 
satisfies an estimate 

Wip{z) = 0(e~Ti zTMz ) for \z\ -» oo 

where z = (x,p) then the symplectic capacity of W : z T Mz < h must be 
>\h. 

Actually, we also state a version of Hardy's theorem, which is valid for 
an arbitrary pair of Lagrangian frames, i.e a transversal pair of Lagrangian 
planes. Therefore, the main results of our investigation provides a rigorous 
justification of a reformulation of the uncertainty principle in [13] due to 
Guillemin and Sternberg: "The smallest subsets of classical phase space in 
which the presence of quantum particle can be detected are its Lagrangian 
submanifolds." Consequently, one could say that one of the main aims of 
the present article is to exploit the symplectic nature of Hardy's uncertainty 
principle in the sense of Guillemin and Sternberg. 

Our work is structured as follows: 

• In Section [2] we prove a multi-dimensional variant of Hardy's theorem, 
as a property of the symplectic spectrum of the matrix ex- 
tracted from the conditions ([3]) where A and B are positive-definite 
symmetric matrices; here F denotes the /i-dependent n-dimensional 
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Fourier transform. In Lemma [T] (Subsection I2.1j) we show that it is 
possible to perform a symplectic diagonalization of a positive-definite 
block-diagonal matrix using symplectic block-diagonal matrices. 

In Section we give a purely geometric interpretation of Hardy's un- 
certainty principle in terms of the notion of the symplectic capacity, 
which is closed related to Gromov's non-squeezing theorem; we take 
the opportunity to quickly review the main definitions and properties 
concerning these objects. In particular, we point out that all symplec- 
tic capacities agree on phase-space ellipsoids; 

In Section 0] we restate the results above in terms of the Wigner trans- 
form, by showing that the conditions © are equivalent to 

WiP{x,p) < Ce~^ xTAx+pTBp) (4) 

for some constant C > 0. In Subsection 14.31 we give an equivalent 
geometric statement of the results above in terms of the topological 
notion of symplectic capacity. We will in fact prove that if ([1]) holds 
for ifi ^ then the symplectic capacity of the phase space ellipsoid 
W : \z t Yj~ 1 z < 1 is at least hh; here S is the covariance matrix 
defined by 

'(Ax) 2 

o (A P y 



where Ax = ^/fijla and Ap = wh/2b. The condition ab < 1 is thus 
equivalent to the Heisenberg inequality AxAp >\h. 

• In Section [5] we give two non-obvious extensions of the results obtained 
in the previous sections. The first extension (Subsection 15. 1|) consists 
in replacing the x,p coordinate system by an arbitrary "Lagrangian 
frame" (£, £') and to use the transitivity of the action of the symplectic 
group on the set of all such frames. In the second extension (Subsection 
I5.2p we consider estimates of the type Wip(z) < Ce - ^^ where Q is 
a twice continuously differentiable function which is uniformly convex. 
We express the necessary condition on that function in terms of the 
symplectic capacity of the convex set Q(z) < h. 

Notation. We will use the shorthand notation z = (x,p) for points of 
the phase space M? n = R n x W 1 . The symplectic product of two vectors 
z = (x,p), z' = (x',p') in R 2n is 

a(z, z) = p ■ x — p ■ x 
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where the dot • stands for the usual (Euclidean) scalar product; alternatively 

a(z, z') = Jz-z' where J = f ^ ^\ (0 (resp. I) is the zero (resp. identity) 

matrix of order n). The symplectic group is denoted by Sp(n): we have 
S G Sp(n) if and only if S is a real matrix of order 2n such that a(Sz, Sz') = 
a(z, z')\ equivalently S T JS = SJS T = J. 

When M is a symmetric matrix we will often write Mx 2 , Mp 2 , Mz 2 
instead of Mx-x (or x T Mx), Mp-p, Mz-z. To express that M is symmetric 
and positive-definite we will use the notation M > 0. 

F denotes the n-dimensional ^-dependent Fourier transform. It is the 
unitary operator L 2 (R n ) — > L 2 (R n ) defined for ip G 5(lR ra ) by 

F ^(P) = i^nY' 2 f e-^(x)d n x. (5) 
2 Hardy's Theorem in Dimension n 

Using classical results on the simultaneous diagonalization of a pair of sym- 
metric matrices it is possible to extend Hardy's theorem to the case of W n 
(see for instance Sitaram et al. [22] ). We are going to prove a variant of 
this result using a symplectic diagonalization; this will allow us to relate our 
statements to the notion of symplectic capacity later on in this work. 



2.1 A symplectic diagonalization result 

The following result, although being of an elementary nature is very use- 
ful. We will see that it is a refined version of Williamson's diagonalization 
theorem [23] in the block-diagonal case. 

We make the preliminary observation that if A and B are positive definite 
matrices then the eigenvalues of AB are real because AB has the same 
eigenvalues as the symmetric matrix A l l 2 BA l l 2 . 

Lemma 1 Let A, B > 0. There exists L G GL(n,M) such that 

L T AL = L~ l B(L T )~~ l = A (6) 

where A = diag(\/Aij ■■■,V^~n) is the diagonal matrix whose eigenvalues are 
the square roots of the eigenvalues Ai,...,A n of AB. 

Proof. We claim that there exists R G GL(n,M) such that 

R T AR = I and R~ 1 B(R T )~ 1 = D (7) 
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where D = diag(Ai, A n ). In fact, first choose P £ GL(n,M) such that 
P T AP = I and set B^ 1 = P T B~ 1 P. Since B^ 1 is symmetric, there exists 
H £ 0(n,R) such that B^ 1 = H T D~ l H where D' 1 is diagonal. Set now 
R = PH T ; we have R T AR = I and also 

R- 1 B(R T )~ 1 = HP- 1 B(P T Y 1 H T = HB\H T = D 

hence the equalities (|7|). Let A = diag( v / Ai, ^^K l )■ Since 

R T AB{R T )- 1 = R T AR(R- 1 B(R T )- 1 ) = D 

the diagonal elements of D are indeed the eigenvalues of AB hence D = A 2 . 
Setting L = RA l > 2 we have 

L T AL = k l l 2 R T ARK 1 / 2 = A 
L- x B{L~ x f = K-^R^BiRTy 1 ^ 1 ' 2 = A 

hence our claim. ■ 

The result above is a precise statement of a classical theorem of Williamson 
|24j in the block-diagonal case. That theorem says that every positive- 
definite symmetric matrix can be diagonalized using symplectic matrices. 
More precisely: let M be a positive definite real 2n x 2n matrix; the eigen- 
values of JM are those of the antisymmetric matrix M 1 / 2 JM 1 / 2 and are 
thus of the type ±«AJ with AJ > 0. We have: 

Theorem 2 (Williamson) (i) There exists S £ Sp(n) such that S 7 MS = 

(o a) where A = diag ( A i' -•' A ™)- 

(ii) The symplectic matrix S is unique up to a unitary factor: if S' is 
another Williamson diagonalizing symplectic matrix then S(S')^ 1 £ U(n). 

Proof, (i) See for instance [6l [IB] for "modern" proofs, (ii) See de Gosson 
[6j- ■ 

We will always arrange the AJ in decreasing order: X" > Af > • • • > A^ 
and call (AJ, A^) the symplectic spectrum of the positive definite matrix 
M. The positive numbers AJ (which only depend on M, and not on S) are 
the Williamson invariants of M. Writing the diagonalizing symplectic ma- 
trix as S = (Xi, ...,X n ; Y 1; Y n ) where the Xj and are column vectors, 
the set B = {Xi, ...,X n ;Yi, ...,Y n } is called a Williamson basis for M (it is 
of course not uniquely defined in general). A Williamson basis is a symplec- 
tic basis of (M. 2n ,a), that is a(Xj,Xj < .) = a(Xj,X/ < .) and a(Yj,Xk) = 5jk for 
1 < j, k < n. 

The following result relates Lemma [1] to Williamson's theorem: 
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Lemma 3 Let A,B > 0. The symplectic spectrum (AJ,...,A^) of M 
A 0\ 

consists of the decreasing sequence > • • • > V%i °f square 



£ y 

roots of the eigenvalues Xj of AB. 

Proof. Let (AJ, A^) be the symplectic spectrum of M. The AJ are the 
eigenvalues of 

they are thus the moduli of the zeroes of the polynomial 

P(t) = det(t 2 / + AB) = det(t 2 / + D) 

where D = diag(Ai, A n ); these zeroes are the numbers ±iyAj, j 
l,...,n; the result follows. ■ 
For L invertible set 

'L~ x 



Ml ={o if)- (8) 

Obviously Ml G Sp(n); Lemma [U can be restated by saying that if (A, B) 
is a pair of symmetric positive definite then there exists L such that 

Lemma[T]is thus a precise version of Williamson's theorem for block-diagonal 
positive matrices — it is not at all obvious from the statement of this the- 
orem that such a matrix can be diagonalized using only a block-diagonal 
symplectic matrix! 



2.2 Application to Hardy's theorem 

Lemma Q] allows us to give a simple proof of a multi-dimensional version of 
this theorem. The following elementary remark will be useful: 

Lemma 4 Letn>l. Fori < j < n let fj be a function of (x±, ..,Xj, ...,x n ) G 
M™" 1 (the tilde~ suppressing the term it covers), and gj a function of Xj G R. 
If 

h = fl 8> gi = ■ ■ ■ = f n ® 9n 

then there exists a constant C such that h = C{g\® ■ ■ ■ <g> g n ). 
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Proof. Assume that n = 2; then 

h(x 1 ,x 2 ) = f 1 (x 2 )gi{x 1 ) = f 2 {xi)g 2 {x 2 ). 

If gi{xx)g 2 (x2) / then 

h(x 2 ) / g 2 (x 2 ) = f 2 {xi) / gi{xi) = C 

hence fi(x 2 ) = Cg 2 (x 2 ) and h(x 1 ,x 2 ) = Cg 1 (xi)g 2 (x 2 ). If gi(xi)g 2 (x 2 ) = 
then h{xi,x 2 ) = hence h(x\,x 2 ) = Cgi(xi)g 2 (x 2 ) in all cases. The general 
case follows by induction on the dimension n: suppose that 

h = fi (8> gi = ■ ■ ■ = f n ® g n = fn+i ® g n +i; 
for fixed x n+ \ the function k = fi (8> g\ = ■ ■ ■ = f n ® g n is given by 

k(x,x n+ i) = C(x n+ i)gi(xi) ■ ■ ■ g n (x n )- 
Since we also have 

k(x,x n+ i) = f n +l(xi, ...,x n )g n+ i(x n+ i) 
it follows that C(x n+ i) = C. ■ 

Theorem 5 Let A and B be two real positive definite matrices and ip G 
L 2 (R n ), tp / 0. Assume that 

\4>(x)\ < C A e~^ Ax2 and \Fip{p)\ < C B e~^ Bp2 (10) 

for some constants Ca,Cb > 0. Then: 

(i) The eigenvalues Xj, j = l,...,n, of AB are < 1; 

(ii) If Xj = 1 for all j, then ip(x) = Ce~2K Ax2 for some some complex 
constant C . 

_J_ Ax 2 

(Hi) If Xj < 1 for some j then ip(x) = Q(x)e 2ft x for some polynomial 
function Q :R n — ► C. 

Proof. Proof of (i). It is of course no restriction to assume that Ca = Cb = 
C. Let L be as in Lemma Q] and order the eigenvalues of AB decreasingly: 
X\ > X 2 > ■ • • > X n . It suffices to show that Ai < 1. Setting iI>l(x) = ^{Lx) 
we have 

f^(p) = f^((l t )- 1 p); 
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in view of (|6|) in Lemma Q] condition (jlOp is equivalent to 

\^l(x)\ < Ce~^h hx2 and \F^ L {p)\ < Ce'^ 2 (11) 
where A = diag(A 1; A2, A n ). Setting tpL,i(xi) = iPl(xi, 0, 0) we have 

\4>L,i(xi)\<Ce~^ x2 . (12) 
On the other hand, by the Fourier inversion formula, 

J F^ L (p)dp 2 ■ ■ ■ dp n = (2irh) n / 2 J j e l i/> L (x)dxdp 2 ■ ■ ■ dp n 
= (2vr/ l )(™- 1 )/ 2 F^,i(Pi) 

and hence 

\F^ L ,i(pi)\ < C Lil e~^ p2 (13) 

for some constant Cl.i > 0. Applying Hardy's theorem to the inequalities 
(fT2~j) and (fl~3j) we must have Af < 1 hence the assertion (i). Proof of (ii). 
The condition Xj = 1 for all j means that 

\iI>l(x)\ < Ce~^n x2 and \Fi/) L (p)\ < Ce~^ p2 (14) 

for some C > 0. Let us keep x' = (x 2 , ■ ■■■,x n ) constant; the partial Fourier 
transform of ipL in the x\ variable is Fiipi = {F')~ 1 FipL where {F')^ 1 is 
the inverse Fourier transform in the x' variables, hence there exists C > 
such that 

n-l f 1 2 

\Fnl> L (x u x')\ < (^) 2 / \F^ L (p)\dp 2 • • • dp n < C'e~2E p K 

12 1/2 

Since \iPl(x)\ < C(x )e~^ Xl with C(x ) < e~^ x it follows from Hardy's 
theorem that we can write 

Mx) = h(x')e^ x2 

for some real C°° function f\ on M n_1 . Applying the same argument to the 
remaining variables x 2 , —,x n we conclude that there exist C°° functions fj 
for j = 2, n, such that 

_ _j_ 2 

4>l(x) = fj(xi, ..,xj, ...,x n )e 2h x i. (15) 
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In view of Lemma U] above we have = Ci,e~^ x for some constant 

Cl; since A = I = L T AL we thus have ip(x) = Cie~ Ax ' 2fi as claimed. 
Proof of (Hi). Assume that Ai < 1 for j G J , J a subset of {1, ...,n}. By 
the same argument as in the proof of part (ii) establishing formula (|15p , we 
infer, using Hardy's theorem in the case ab < 1, that 

^l(x) = fj(x 1 ,..,x j ,...,x n )Q j (x j )e ih x i 

where Qj is a polynomial with degree if j ^ J . One concludes the proof 
using one again Lemma |U ■ 



3 Geometric interpretation 

Let us give a geometric interpretation of Theorem [5j We begin by making 
an obvious observation: Hardy's uncertainty principle can be restated by 

1 2 

saying that if ip ^ then the conditions ip(x) = 0(e~^ ax ) and Fijj(p) = 
0(e~2h bp ) imply that the ellipse W : ax 2 + bp 2 < h has area nh/Vab > 
irh =\h\ 

Area(W) >-h. 



More precisely: 

// the area of the ellipse W is smaller than \h then ip = 0; if this area 
equals \h then i\)(x) = Ce~^ ax2 and if it is larger than \h then ip(x) = 

1 2 

Q(x)e~2ft ax where Q is a polynomial function. 

When trying to generalize this observation to higher dimensions, one 
should resist the pitfall of copying the statement above mutatis mutandis 
and replacing everywhere the word "area" by "volume". As we will see, 
volume is not the right answer; one has instead to use the more subtle 
notion of symplectic capacity, introduced by Ekeland and Hofer [3j following 
Gromov's work [13] on pseudoholomorphic curves. 



3.1 Symplectic capacities and symplectic spectrum 

A symplectic capacity on the symplectic space (M 2n , a) assigns to every 
subset Q of M 2n a number c(O) > or +oo; this assignment has the four 
properties listed below. (We denote by B(R) the ball \z\ < R and by Zj(R) 
the cylinder x 2 + p 2 < R 2 .) 

SCI Monotonicity: c(VL) < c(fi') if n C n'; 
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SC2 Symplectic invariance: c(/(0)) = c(f2) for every symplectomor- 
phism f defined near f2; 

SC3 Conformality: c(Afi) = A 2 c(ft) if X £ R; 

SC4 Nontriviality: We /iaue c(B(R)) = c{Zj(R)) = ttR 2 . 

A fundamental example of symplectic capacity is provided by the "Gro- 
mov width" , defined by 

c Gr (Q) = sup {irr 2 : f(B(R)) C $7} (16) 

/eSymp(n) 

where Symp(n) is the group of all symplectomorphisms of (M. 2n ,a). Prop- 
erties (i)-(iii) and cq t (B(R)) = irR 2 are trivially verified; that we also have 
CQ r {Zj{R)) = irR 2 is just Gromov's non-squeezing theorem [13] which as- 
serts that a phase-space ball cannot be squeezed inside a symplectic cylinder 
with smaller radius using symplectomorphisms (but such a squeezing can, 
of course, be performed using general volume-preserving diffeomorphisms) . 
We will also use the "linear symplectic capacity" c nn defined by 

cun(n)= sup {nr 2 : f(B(R)) C Q} (17) 

/elSp(n) 

where / this time ranges over the group ISp(n) of all affine symplectic au- 
tomorphisms of (M 2n , a) (the "inhomogeneous symplectic group" ) . The ca- 
pacity cii n has the same properties as general symplectic capacities, except 
that it is only invariant under linear or affine symplectomorphisms. 

We have the following result, which allows us to talk about the symplectic 
capacity of a phase-space ellipsoid: 

Lemma 6 For M > let Q M = {z € R 2n : Mz 2 < 1}. For any symplectic 
capacity c on (M? n ,a) we have 

c(flM) = CKn (fi M ) = ^ (18) 
A i 

where A J > • • • > is the symplectic spectrum of M. 

Proof. See for instance Hofer and Zehnder [16], Proposition 2, §2.1, p. 54 
or de Gosson [6], Proposition 8.25, p. 251 (where A^ should be replaced by 
XI). m 
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3.2 Application to the Wigner ellipsoid 

We can restate Hardy's theorem in a very simple geometric way in terms of 
the symplectic capacity of the "Wigner ellipsoid" (the terminology seems to 
be due to Littlejohn |20j): 

Proposition 7 Let tp 6 L 2 (R n ), if) / 0. Assume that 

\^{x)\<C A e-^h Ax2 and \F^(p)\<C B e-^i Bp2 . (19) 
Then the symplectic capacity of the Wigner ellipsoid 

W : Ax 2 + Bp 2 < h 

satisfies c(W) > \h. 

Proof. Setting M = the equation of W is Mz 2 < h. Consider the 

ellipsoid Om : Mz 2 < 1. Let (Aj , Af, A^) be the symplectic spectrum of 
M; by formula (|18p in Lemma[6]we have c(£Im) = ^/X^. In view of Lemma 
|3] AJ = \f~X~j where the Xj are the eigenvalues of AB, and by Theorem 
[5] we must have Xj < 1, hence c(f2j\f) > tt- Since we have 

c(W) = Hc(Q,m) in view of the conformality property (CZ3); the result 
follows. ■ 

The result above will be extended to the Wigner distribution in next 
section. 

4 Hardy's Theorem and Wigner's distribution 

It turns out that Hardy's theorem - which involves two conditions, one 
about a function and the other about the Fourier transform of that function 
- is equivalent to a single condition on the Wigner transform of ip. This 
condition will be made explicit in Theorem [9] below; let us first prove some 
preliminary results about Wigner transforms. 

4.1 Wigner Distributions 

The Wigner transform of a function was introduced by Wigner in |23j . fol- 
lowing joint work with Szilard. It is defined, for ip £ L 2 (R n ), by the formula 

W^iz) = (^) n J e- l n^{x + \y)^{x-\y)<Fy. (20) 
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The function Wip has a simple interpretation in terms on the theory of Weyl 
pseudodifferential operators. For ip £ L 2 (M n ) with = 1 consider the 
orthogonal projection P^ on the ray {A^ : A £ C}; we have P^<p(x) = (<f>\ip)ip 
for <p £ L 2 (M. n ) hence the operator kernel of P^ is = tp <8>ip. Writing P^ 
in Weyl operator form we have 

P^(x) = jj e^< x -^p^\{x + y),p)ct>{y)d n pd n y 

where the symbol p^ is given by 

Pf(x,p) = (^h)" K<p(x + \y,x- ~y) = Wtp(z). 

More generally, one might want to consider the cross- Wigner transform 
(also called Wigner-Moyal transform) which associates to a pair (ip,4>) £ 
L 2 (R n ) x L 2 (M n ) the function 

W^,cf>)(z) = I e-^{x + \y)4>{x-\y)^y; 

(it is the Weyl symbol of the operator defined by the kernel ip®4>); of course 
W(ijj,(j)) = Wip. The following properties of the (cross-) Wigner transform 
are well-known 

Wl W(if), <f>) = W(<t>,ijj) (hence Wip is real); 
W2 If if;, Ftp £ L 1 (M n ) n L 2 (R n ) then 

I W*p(z)d n p = |^(x)| 2 , I W*p(z)d n x = \F^(p)\ 2 . (21) 

Recall that the metaplectic group Mp(n) is generated by the following 
unitary operators on L 2 (R 2n ) (OS]): the scaling operators ML jTn (L £ 
GL(n,R)), the "chirps" Vp {P = P T ), and the modified Fourier transform 
J = i~ n l 2 F\ by definition 

M L ^(x)=i m ^\^L\^{Lx) , V P ^(x) = e^ Px2 iP(x) (22) 

(m corresponds to a choice of argument for detL). Mp(n) is a faithful 
representation of the double covering group of Sp(n); the projection tt : 
Mp(n) — * Sp(n) is determined by its action on the generators: 

vr(M Lim ) = M L , tt(Vp) = (_ J p fj , tt( J) = J (23) 

(M L defined by formula ©). 
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W3 Let S be any of the two metaplectic operators associated with S £ 
Sp(n). The following metaplectic covariance formula holds: 

W(Sip)(z) = Wip^z). (24) 



As a particular case of (|24|) we have 

W(Fijj)(z) = W{J^)(z) = Wip(-Jz). (25) 

Recall that the Heisenberg operator T{zq) is defined, for zq = (xq,pq) £ 
R 2n , by 

f(zoMx) = e^ p °- x -^°- Xo) ij(x - x ). (26) 

W4 We have 

W(T(z )iP)(z)=Wi>(z-z ). (27) 

The Wigner transform behaves well under tensor products: if x = {x\ x") 
with x' £ R k , x" £ M. n ~ k and ip 1 £ L 2 (R k ), ip" £ L 2 (R n - fc ), then 

W(Tp' (8 ip") = W'tp' W"ij}" (28) 

where W and W" are the Wigner transforms on L 2 (R. k ) and L 2 (W n ~ k ), 
respectively. More generally, if W,W',W now denote cross- Wigner distri- 
butions: 

W(ip' ® V", 8> ^") = W'(ip', <f/) <g> 0")- (29) 



4.2 Wigner transform and Hermite functions 

The A:-th state of the quantum harmonic oscillator with classical Hamiltonian 
H(x,p) = \{x 2 + p 2 ) is the Hermite function 

_j_ 2 

{hk the fe-th Hermite polynomial). One shows that 

rmn(k,£) 

W(^ k ,^)(z)=e-Ti^ 2 Yl C.ikJ)/-^ (30) 

where the Cj(k,£) are real constants and z is identified with x + ip £ C n in 
the right-hand side (see e.g. [4], p. 66-67). Notice that in particular 

\Wty k ,4> e ){z)\ <e- l h\ z \ 2 P M (\z\) (31) 
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where Pki is a real polynomial of degree k + £. 

We will need the following Lemma which says that the Wigner transform 
of a Hermite function is the product of an exponential and of a polynomial 
with positive leading coefficient. (For related results see [18J). 

Lemma 8 Let Q be a (complex) polynomial function on M. n and ip(x) = 
Q(x)e-^ Ax \ A > 0. Then: 

(i) The Wigner transform of tp is given by 



2 



Wip(x,p) = R{z A ,zX)e~ h\ z *\- (32) 

where R is a polynomial function and za = A 1 l 2 x + iA~ l l 2 p (A 1 / 2 the 
positive square root of A); 
(ii) In particular 

\Wip(x,p)\ < T(\z A \)e-h z *\ 2 (33) 
where T is a polynomial with real coefficients. 

Proof, (i) Let us set cp = M A -i/2 ip where M A -i/ 2 £ Mp(n) is defined by 
(|22jh Thus 

ip{x) = P{x)e-^r>\ 2 with P{x) = V det A~ 1 Q(A' 1/2 x) 



and we have, by property (|24l) of the Wigner transform and the first formula 

dSSD, 

Wij(z) = W(p(A 1/2 x, A~ 1/2 p). (34) 

Writing P(x) = a a x a (we are using multi-index notation a = (a±, ...,a n ), 
x a = x" 1 • • • x° n ) we have 

<^(x) = V a a ip a (x) , <p a = vV®---®Vn n 

with ip^ (xj) = x^ e~ x ^ 2 . By the sesquilinearity of the cross- Wigner trans- 
form we get 

W<p = J2 a a a-pW{ip a ,^) (35) 
and by the tensor product property (|29p 



W(<p a , <fP) = W(<p? ® ■ ■ • ® <p%>, ¥>f ® ■ ■ • ® <pfr) 
= W(^ 1 ,^ 1 )8)---(8)W«» ) ^«). 
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The Hermite functions tpk forming an orthonormal basis of L 2 (R) each yy 7 
is a finite linear combination of these functions; using again sesquilinearity 
and applying formula (|30p there exist polynomials Paj,fy such that 

with zj = Xj + ipj and hence 

W(^,^)(z)=P aP (z,z)e~h^ 

where P a p = P ai ,pi ® • • • <8> -Pa„,/3 n is a polynomial function in 2n variables. 
It follows from ([35 p that 

W^(ar) = ^ a/3 a a o^P Q/3 (z,z)e-ll 2 l 2 = i?(z,z)e-ll z l 2 (36) 

and hence, in view of ([341) . 

= E „ w^(^ 1/2 ^> A- 1 / a p) e -i^- 1 * a +^) 

as claimed, (ii) Since T^y? is a real function we have Wip(z) < \Wip\ and 
hence, taking (|36p into account, 

WV(z) < |i?(z,^)|e-Hl z l 2 < T(\z\)e~h z \ 2 
which concludes the proof in view of (|34|) . ■ 

4.3 Phase-space formulation of Hardy's theorem 

When dealing with Gaussian functions related to "squeezed coherent states" 
we obtain Gaussian estimates where the quadratic form in the exponent no 
longer is block-diagonal. For instance, the Wigner transform of a Gaussian 
of the type 

tyx,Y\x) = e 2h y ' 
(X and Y real symmetric, X > 0) is given by the formula 

W^ x ,y{z) = {Kh)- n l 2 {&etX)- 1 l 2 e-^h Gz2 (37) 

where the matrix G is given by 

(X + YX~ l Y YX~ 
G ~ V X- l Y X- 1 




16 



(see Proposition 8.4, p. 263, in de Gosson [6]); the result seems to go back 
to Bastiaans according to Littlejohn [20]). An important observation is 
that G is a positive-definite symplectic matrix as follows from the obvious 
factorization 

G = S T S wtth S=L-i/3y x-V*) € Sp(n) " (39) 

Setting S = \G- X the ellipsoid W : |5TV < 1 is the set {z : S T Sz 2 < h}; 
W is thus the image of the ball B(yK) by a linear symplectic transformation, 
and thus has symplectic capacity \h. 

Let us now show, as claimed in the introduction, that Hardy's uncer- 
tainty principle for a function ip is equivalent to a condition on its Wigner 
transform Wip. 

Theorem 9 Let ip S L 2 (W l ) and A, B two positive real nxn matrices. Let 
Ca,Cb > 0. The condition 

\ip(x)\ < C A e~^ Ax2 and \Fi/>(p)\ < C B e~^ Bp2 (40) 
is equivalent to the existence of a constant Cab > such that 

Wi){z) < C AB e-^ {Ax2+Bp2) . (41) 

Proof. In view of properties (|2ip of Wip, condition (I41D implies that there 
exist constants Ca, Cb > such that 

hK*)| 2 < tikrW , \F^{p)\ 2 <C 2 B e~Ti Bp2 

hence ([4T ]) ^> pD]) . Let us prove that conversely (@0])=>(|II]). Let Ai, \ n 
be the eigenvalues of AB. If there exists j such that Xj > 1 then ip = 
by Theorem [5] and (|4ip is trivially verified. We may thus assume from 
now on that Xj < 1 for j = 1, ...,n. Let L S GL(n,M) be as in Lemma [H 
that is L T AL = L~ 1 B(L T )~ 1 = A where A is the diagonal matrix whose 
eigenvalues are the y/Xj. We have, setting iPl(x) = ip(Lx) as in the proof 
of Theorem [5j 

\i>Lix)\ < C A e~^ Kx2 and \F^ L (p)\ < C B e-^ Ap2 . 
Since Xj < 1 for all j = 1, n Theorem [5] implies that we have 

iPl(x) = QL{x)e~h hx2 
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where Ql is a polynomial function which is constant when all the Xj are 
equal to one. It follows, by Lemma [8] that 

WM*) < R L (yi-,yn)e-^ Ax2+A ~ lp2) 

where Rl is a polynomial function with positive leading coefficient and yj = 
XjX 2 + XJ p 2 . Let Cl > be a constant such that 

n , 

for every e > there exists Cl, £ > such that 

n 

J]^' e -l( A - 2 + A -V) < Ci)£e -i((A- £ )^+A-V) 

i=i 

(we are writing A — e for A — Is) and hence 

W^ L {z) < C L , £ e-H(( A -^' 2 + A "V). (42) 
Applying the same argument to W(FipL)(x,p) = WipL(—p,x) we also have 

W^ L (z) < C Lt£ e-TS Ax2 +( A ' 1 ~^ 2 \ (43) 

Since 

sup[(A - e)x 2 + A~ 1 p 2 ,Ax 2 + (A -1 - e)p 2 } = Ax 2 + A~V 
the inequalities (j4"2j) - (j4'3"|) imply that we have 

WMz)<C L ,se-^ Ax2+A ~ lp2) -, 

since iPl(%) = 4>{Lx) this is just condition (|4T|) . ■ 

Theorem [9] has the following consequence which contains Hardy's the- 
orem as a particular case (we have proved a particular case of that result, 
using different methods, in [U [9]). 

Corollary 10 Let tp £ L 2 (W l ), ^ 7^ 0. Assume that there exists a positive- 
definite real matrix M, a vector a E M? n and C > such that 

Wi>{z) < Ce -Ti {Mz2+2a - z \ (44) 

Then the ellipsoid W = {z : Mz 2 < h} has symplectic capacity c(W) > ^h 
(equivalently X " < h). 
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Proof. Assume first that a = 0. Let S G Sp(n) be such that S T MS is in 
Williamson diagonal form 

'A N 



° -0 A 

with A = diag(Af, A£), Xf > ••• > A£. Choose 5 G Mp(n) with projection 
S. It follows from the metaplectic covariance property (|24|) of the Wigner 
transform that 

Applying Theorem ([9]) there exist constants C\ , C2 > such that 

\S~ l ip(x)\ < C ie -^ Ax2 and \FS- l ^{x)\<C 2 e'^h Ap2 . 

In view of the multidimensional Hardy theorem we must have AJ < 1 for 
j = 1, ...re hence (Lemma|6| c(W) = irh/\1 > 7r?t which concludes the proof 
in the case a = 0. Assume now a is arbitrary, and set Q(z) = Mz 2 + 2a • z; 
choosing S G Sp(n) and 5 G Mp(n) as above we have Q{Sz) = Dz 2 + 2b ■ z 
where b = S T a; completing squares we get 

Q(Sz) = D(z + D~ x bf - D- l b 2 . 

It follows that for a new constant C we have 

W(S-^)(z) < Ce-l Q ^ < C'A D ^ D ~^ 2 . 

We next observe that 

W(S~ l i)){z - D^b) = W(T(b)S- x il))(z) 

where T(D~ 1 b) is a Heisenberg operator; we thus have, using ([2 



W{T{D^ 1 b)S^ lr 4))(z) < C'e-7i D 



and it now suffices to apply the case a = to ip' = T(D 1 b)S tp. ■ 

It is instructive to see how the sub-Gaussian estimate (1441) is related to 



the uncertainty principle of Quantum Mechanics. Setting £ = 1 we 
can define the multivariate Gaussian probability density 

P (,) = (i-) n (detS)- 1 /V^ 2 

and view £ as a statistical covariance matrix. The Wigner ellipsoid W : 
Mz 2 <his identical with the set W s = {z : §£ _1 z 2 < 1}. 
One of us has proven in [B] the following result: 
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Proposition 11 The two following conditions are equivalent: 

(i) c(W) = c(W s ) > \h 

(ii) The Hermitian matrix E + ii J is positive semidefinite. 



Write now E in block-matrix form 

E 



f^xx Exp 
\Epx Epp 



where S^i, E^p = Ep^, an d Epp are the n x n partial covariance ma- 
trices Sxx = (Cov(a;j,Xfc))j i fc, T, X p = E^ P = (Cov^pfc))^, and Epp = 
(Cov (pj,pk))j,k] the covariances are defined with respect to the probability 
density : setting Zj = Xj for 1 < j < n and Zj = pj for n + 1 < j < 2n we 
have 

Cov(zj,Zfc)= y ZjZj t p{z)d 2n z — j Zjp(z)d 2n z J z^p{z)d 2n z 

The conditions (i) and (ii) in Proposition 1111 above are equivalent to 

(A Xj ) 2 (A P:j ) 2 > (Covix^pj)) 2 + \h 2 (45) 

for j = l,...,n (see Narcowich [21] and de Gosson [6], and the references 
therein) . The inequalities (|45p (known in the quantum-mechanical literature 
as the Schrodinger-Robertson uncertainty relations) are a precise form of 
the usual text-book Heisenberg inequalities AXjAPj > to which they 
reduce if one neglects correlations. 



5 Two Extensions of Hardy's Theorem 

5.1 Restatement in an arbitrary Lagrangian frame 

The statements of Hardy's uncertainty principle we have been considering 
correspond to a particular choice of coordinates namely the positions x and 
the momenta p for the phase space. These statements thus correspond to 
the choice of frame (£ x ,£p) where tx is the horizontal Lagrangian plane 
R n x {0} and lp the vertical Lagrangian plane {0} x M. n . This choice is of 
course to a great extent arbitrary. In the following we are going to extend 
our results to arbitrary Lagrangian frames. 

Recall that a subspace £ of the phase space (R 2n , a) is called isotropic, if 
the symplectic form a vanishes identically on £. If £ has maximal dimension 
n, then £ is called a Lagrangian plane. The set of all Lagrangian planes in 
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(M? n ,a) is called the Lagrangian Grassmannian of (IR 2 ™, a) and denoted by 
Lag(n). 

The subgroup of all symplectic matrices S such that St = £ is called 
the stabilizer of £ and denoted by St (I). Note that S G St(£) if and only if 



If £ and £' are Lagrangian planes in (M 2n , a) satisfying £ n £' = 0, then 
are called transversal; equivalently £(&£' = M. 2n . We will call a pair 
of transversal Lagrangian planes a Lagrangian frame. An important 
property is that the symplectic group Sp(n) acts transitively not only on the 
Lagrangian Grassmannian Lag(n), but also on the set of Lagrangian frames: 
if (£i,£i) and (£2, £'2) are pairs of Lagrangian planes satisfying £\ n £\ = 
£ 2 n£' 2 = 0, then there exists S £ Sp(n) such that (£ 2 ,£ 2 ) = (S£i,S£[) (see 
de Gosson [6]). 

We interpret the marginal properties (|2ip of Wigner's distribution in 
terms of the horizontal and vertical Lagrangian plane lx and £p: if t/j, Ftjj £ 
L 1 (R n ) n L 2 (R n ) then we can rewrite (HU) as 



Recall that the metaplectic covariance property (W5) of the Wigner distribu- 
tion tells us that if S £ Mp(n) has projection S on Sp(n) then W(Sip)(z) = 
Wip(S~ 1 z). Therefore, we get 



Since Jlx = £p and J = F we get the analogous results for the Lagrangian 
plane SJ. Note that {S£x , SJix) are a transversal pair of Lagrangian 
planes. In other words the transitivity of the symplectic group on Lag(n) 
allows to translate a statement about £x and £p into a statement of another 
pair of Lagrangian planes (£,£')■ One just has to choose the correct S to go 
from (£x,£p) to (£,£')■ Then the statements about tp and Ftp translate into 
statements about Sip and S o Jqp. Consequently, one of our main results, 
Theorem [9l remains valid in an arbitrary Lagrangian frame, if one makes 
the proper modifications as indicate above. 



S T g St(J£). 





or, equivalently, 
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5.2 The case of convex exponents 

We are going to extend the Corollary [TU] of Theorem to the case where the 
inequality (|4"4"|) is replaced by 

Wip(z) < Ce'^ Q{z) (46) 

where Q is a uniformly convex function on R 2n ; we will assume that Q(0) =0 
(the case Q(0) ^ is trivially reduced to this case by changing the constant 
C). Using the same trick as in the proof of Corollary 1101 we may moreover 
assume, replacing by T(zo)V' f° r a suitably chosen zq £ M. 2n , that 

Q'(0) = V z Q(0) = 0. 

Let us briefly recall a few basic facts on convex functions (see Andrei 
PQ for a concise review of the topic). A function Q : R 2n — > 1R is strictly 
convex if we have 

Q(az + (1 - a)z) < aQ(z) + (1 - a)Q(z') 

for < a < 1 and z ^ z' . If the function Q is of class C 2 this condition 
is equivalent to Q"(z) > for all z S M 2n (Q"(z) is the Hessian matrix of 
Q calculated at z). In what follows we will make in addition the following 
uniformity assumption: 

There exists c > such that Q"{zq)z 2 > c\z\ 2 for all z,zq 6 M. 2n . (47) 

Diagonalizing Q" {zq) by an orthogonal matrix H{zq) this condition is easily 
seen to be equivalent to the condition \q > where 

Aq = inf {X(z) : X(z) is an eigenvalue of Q"(z)}; (48) 

z eR 2n 

the smallest constant c for which (|47[) holds is then precisely Aq. 

Proposition 12 Under the same assumptions on Q as above letip G L 2 (M. 2n ), 

ip be such that Wtp(z) < Ce~K®^ for some C > 0. Then the convex 
set 

C = {z£ R 2n : Q{z) < h} (49) 
satisfies c(C) > \h for every symplectic capacity c on (M. 2n ,o~). 
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Proof. Let us begin by showing that we have < Xq < 2. Since the 
uniformity assumption is equivalent to Xq > it suffices to show that Xq < 
2. In view of the mean value theorem we have 

Q(z) = Q(0) + Q'(0) • z + \Q\z')z 2 = \Q"{z')z 2 (50) 

where z' lies on the line segment joining to z\ we have Q"(z')z 2 > Aq|z| 2 
hence 

Q(z) > ±X Q \z\ 2 (51) 

so that 

Wip(z) < Ce~Ti Q{z) < C'e-h^W. 

The symplectic spectrum of XqI consists of the point Xq hence we must 
have Xq < 2 choosing M = \XqI in Corollary 1101 The proposition follows: 
in view of (f5"T|) the condition Q(z) < h~ implies ^Aq|z| 2 < h hence the set C 
contains the ball \z\ 2 < 2H/Xq < H; one concluded using the monotonicity 
of symplectic capacities (Property (SCI)). ■ 

Let C be a compact and convex set. We recall (John [19]) that there 
exists a unique ellipsoid W contained in C having maximal volume. This 
ellipsoid, called the "John ellipsoid" HSJ has the property that 

W C C C z + 2n(W - zo) (52) 

where zq is the center of W. The result above has the following immediate 
consequence: 

Corollary 13 Let W be the John ellipsoid associated to the convex and 
compact set ( f^P| ). We have c(W) > \h. 

Proof. The uniform convexity of Q implies that the convex set C = {z : 
Q(z) < h} is compact (Andrei pQ); John's ellipsoid is thus well-defined. In 
the proof of Proposition [T2] we have seen that C contains the ball \z\ 2 < 
2H/Xq < H; this ball is contained in John's ellipsoid. The result follows 
again in view of the monotonicity of a symplectic capacity. ■ 

Proposition [12] has another interesting non-trivial consequence. In [16] 
Hofer and Zehnder construct a symplectic capacity chz having the following 
property: 

If Q is a convex and compact subset of R 2n with smooth 
boundary d£l then 



chz(^) = inf I <b pda 
T U 7 
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where 7 ranges over the set of all periodic Hamiltonian orbits on 
dQ. 



Corollary 14 Under the same assumptions on Q and ifi as above we have 

pdx 



for every periodic Hamiltonian orbit 7 on the hypersurface defined by Q(z) = 
H. 

Proof. The boundary of C is precisely the hypersurface defined by Q(z) = h. 
In view of Proposition [12] we have c(C) > \h for every symplectic capacity 
c hence, choosing c = chz 



11. ■ 1 - 1 



which proves ([53 



inf < <f> pdx > > -h 
1 U7 J 2 



6 Concluding Remarks 

We have proved a n-dimensional version of Hardy's uncertainty principle, 
and showed that it is equivalent to a statement on the Wigner distribution of 
a sub-Gaussian state. The extension of this result to more general estimates 
involving convex exponents in Subsection 15.21 opens the door to the study 
of non-trivial properties for the density matrix of quantum systems. Such 
applications are very important for the understanding of non-linear quantum 
optics and the theory of entangled quantum states. 

We mention that Hogan and Lakey [17] have done a very interesting 
analysis of the interplay between Hardy's uncertainty principle and rota- 
tions. It would certainly be useful to restate their results in our context; we 
leave this possibility for further work. Also, Grochenig and Zimmermann 
|10| have studied Gaussian estimates from the point of view of the short-time 
Fourier transforms; the methods they use are very different from ours. 
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